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The dust scattering model can not explain the shallow X-ray decay 
in GRB afterglows 
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ABSTRACT 

A dust scattering model was recently proposed to explain the shallow X-ray decay (plateau) 
observed prevalently in Gamma-Ray Burst (GRB) early afterglows. In this model the plateau 
is the scattered prompt X-ray emission by the dust located close (about 10 to a few hundred 
pc) to the GRB site. In this paper we carefully investigate the model and find that the scattered 
emission undergoes strong spectral softening with time, due to the model's essential ingre- 
dient that harder X-ray photons have smaller scattering angle thus arrive earlier, while softer 
photons suffer larger angle scattering and arrive later. The model predicts a significant change, 
i.e., A/3 ~ 2 — 3, in the X-ray spectral index from the beginning of the plateau toward the end 
of the plateau, while the observed data shows close to zero softening during the plateau and 
the plateau-to-normal transition phase. The scattering model predicts a big difference between 
the harder X-ray light curve and the softer X-ray light curve, i.e., the plateau in harder X-rays 
ends much earlier than in softer X-rays. This feature is not seen in the data. The large scatter- 
ing optical depths of the dust required by the model imply strong extinction in optical, Ay ;> 
10, which contradicts current findings of Ay — 0.1 — 0.7 from optical and X-ray afterglow 
observations. We conclude that the dust scattering model can not explain the X-ray plateaus. 
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!_< \ 1 INTRODUCTION 



Swift has discovered a generic behaviour in X-ray afterglows of 
Gamma-Ray Bursts (GRB): the X-ray light curve (LC) first shows 
a steep decline during a few hundred seconds after the end of the 
7-rays, then it shows a shallow decay lasting 10 4 — 10 5 s which 
is followed by a "normal" power-law decay (Nousek et al. 2006; 
O'Brien et al. 2006). The normal decay at late times is the canonical 
afterglow component due to the interaction of the decelerated GRB 
ejecta with the circumburst medium, i.e., the forward shock model. 
The steep decline is generally interpreted to have the same origin as 
the prompt 7-ray emission (e.g., Kumar & Panaitescu 2000; Liang 
et al. 2006). 

The intervening shallow decay, sometimes called the 
'plateau', is the most puzzling feature of the X-ray LC. The most 
straightforward interpretation is a late steady energy injection into 
the external shock, where the latter is produced by the decelerated 
early ejecta plunging into the medium. The late energy injection 
could be due to a new ejecta from the late activity of the central en- 
gine (e.g., Dai & Lu 1998a,b; Zhang & Meszaros 2001; Dai 2004; 
Yu & Dai 2007), or due to a slow trailing part of the outflow catch- 
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ing up with the already forward-shock-decelerated early part of the 
outflow when the outflow has a spread in its Lorentz factor distribu- 
tion (e.g., Granot & Kumar 2006). If it is the first scenario, then this 
interpretation implies a steady, late activity of the central engine - 
lasting as long as a day - which poses a challenge to the models of 
the central engine. Moreover, according to the energy-injection in- 
terpretation, the plateau-to-normal transition in the LC corresponds 
to the cessation of the energy injection, thus the transition should 
be achromatic. But in about 1/3 of the X-ray plateau GRBs with op- 
tical afterglow observations, the optical LC does not show a simul- 
taneous plateau-to-normal break, while in another smaller fraction 
of the plateau cases, the plateau-to-normal breaks in optical and 
X-ray are indeed simultaneous (Panaitescu 2007). In most cases 
the power-law decay following the plateau is consistent with the 
predictions (the closure relationships) of the forward shock model, 
which in turn is consistent with the energy injection interpretation. 
There is a long list of alternative models for the plateau phase, such 
as a slow energy transfer from the ejecta to the ambient medium 
(Kobayashi & Zhang 2007), a two-component jet model (e.g., Gra- 
not et al. 2006), a varying shock microphysical parameter model 
(e.g., Panaitescu et al. 2006), and a reverse shock dominated after- 
glow model (Uhm & Beloborodov 2007; Genet et al. 2007), etc. 
(see Zhang 2007 for a review), but none of them satisfy all the ob- 
servational constraints. 

An attractive possibility was suggested by Shao & Dai (2007) 
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regarding the origin of the X-ray plateau. If the long-duration GRB 
progenitors are massive stars, it is very likely that dust exists in the 
vicinity of the GRB site since it is in a star forming region. The 
X-ray photons from the GRB and its afterglow can be scattered in 
small angles by the dust near the line-of-sight to the GRB, as anal- 
ogous to the halo emissions of other X-ray sources (e.g., Smith & 
Dwek 1998). The GRB prompt emission scattered off the dust has 
been considered earlier by Esin & Blandford (2000) and Meszaros 
& Gruzinov (2000). Aside from the scattering by the dust local 
to the GRB site, Miralda-Escude (1999) considered the scattering 
of the X-rays from the GRB afterglows by the dust in the inter- 
vening galaxies along the line-of-sight to the GRB, but the flux 
turns out to be very low and difficult to detect for that case. De- 
pending on the distance of the local dust region to the GRB site, 
a delayed emission component from the scattering can show up in 
the afterglows. Shao & Dai (2007) and Shao et al. (2008) recently 
used this scenario to interpret the plateau phase in the X-ray after- 
glow LC as to be the scattered prompt X-rays by the dust located 
at about ten to a few hundred pc from the GRB site. The scatter- 
ing happens preferentially within a characteristic scattering angle 
6 C which is dependent on the photon energy E and the dust grain 
size. At larger angles the differential scattering cross section of the 
dust grains decays steeply. Therefore the scattering within 9 C gives 
rise to a plateau phase whose duration is determined by 9 C and the 
distance of the dust region to the GRB site. Larger angle scattering 
produces a F(t) oc t~ 2 decay following the plateau. This model 
does not need to invoke a long steady central engine activity. In ad- 
dition, since the scattering only works in the X-ray band, the lack 
of a simultaneous break in optical LC does not pose a problem for 
this model. 

The purpose of this work is to carefully investigate the output 
of this dust scattering model - in terms of the spectral and tempo- 
ral properties of the scattered emission - and to compare it with 
the data. The paper is structured as follows. We first calculate and 
quantify the softening expected from the dust scattering model in 
Section 2. Then, we search in the data for evidence in favour of 
the model including the spectral evolution in the plateau and post- 
plateau phases for a sample of GRBs in Section 3 and 4. An ex- 
pected difference in hard X-ray and soft X-ray LCs is discussed 
in Section 5. We calculate and discuss the optical extinction for 
the dust in Section 6. Our conclusion and further discussion are 
presented in Section 7. Throughout the paper the spectral index j3 
and the time decay index a of the emission flux are defined as in 

/„(t)Ki/-"r. 



2 SPECTRAL SOFTENING IN DUST ECHO EMISSION 

We first derive the temporal and spectral properties of the scattered 
emission or the 'echo' (hereafter we use 'echo' and 'scattered emis- 
sion' interchangeably) by the dust in the simplest geometry where 
the dust is concentrated in a thin layer (or a dust "screen") near 
the GRB, following Shao & Dai (2007). Then we consider a gener- 
alised geometry where the dust is distributed in an extended zone. 

Let us consider a dust "screen" located at a distance R from 
the GRB source. This dust screen does not have to enclose entirely 
the GRB source, as long as its angular size with respect to the GRB 
site is larger than the characteristic scattering angle 8 C (see below). 
The grains in the dust have a size distribution dN(a)/da oc a~ q 
within a range (a_ , a+), where a is the grain size, q is the distribu- 
tion index and N(a) is the column density of all grains with size 
< a. In this paper we use these typical values a_= 0.025 /im, a + = 



0.25 /im, and q= 3.5 inferred from the observations (Mathis et al. 
1977; Mauche & Gorenstein 1986; Draine 2003). We found that 
adopting other typical values did not change our main results. Con- 
sider a GRB source with a fluence per unit energy S(E) [erg cm -2 
keV -1 ] at X-ray photon energy E. Since the GRB source duration 
(~ 10 s) is much shorter than the plateau, it can be considered as 
being instantaneous. 

The flux of the dust scattered emission per photon energy, per 
grain size, at time t can be estimated by 



F B , a (t) 



S(E) 
t ' 
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(1) 



where t[E, a, 6(t)] is the scattering optical depth per grain size a, 
to the photon with energy E and at the scattering angle 9(t); 9(t) 
is given by the geometrical relation t = R8 2 / (2c). 

The angular part of the optical depth can be separated out from 

r by 



-\E,a,0{t)] = 2r a (E)j 2 1 [x(E,a,t)], 



(2) 



where T a (E) is the total optical depth per grain size a and to the 
photon energy E; ji(x) — sin(a;)/a; 2 — cos(a;)/a; is the spherical 
Bessel function of the first order which describes the scattering- 
angle dependence of the cross section, and x = 2irEa9/(hc) is 
the scaled scattering angle where h is the plank constant and c is 
the light speed (Overbeck 1965; Alcock & Hatchett 1978). Via the 
geometrical relation, x can be expressed in terms of E, a and t: 
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ji(x) increases as oc x 2 from x = to x ~ 1.5 and then 
drops rapidly as oc x~ 2 for x > 1.5. Therefore, at a given photon 
energy E, the echo flux LC first appears as a plateau, then tran- 
sitions to a decay as steep as oc t~ 2 . The transition time, which 
corresponds to a characteristic scattering angle 9 C and in turn to 
x ~ 1.5, would be given by 
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We see from Eq. (4) that the duration of plateau is very sensi- 
tive to the photon energy: the plateau at higher energies ends much 
earlier than that at lower energies. Thus the overall echo emission 
must experience strong spectral softening. Note that if the echo is 
observed within a finite energy range, such as in the XRT band (0.3 
- 10 keV), the softening must have begun long before the end of 
the plateau, because the overall plateau ending time is determined 
by t c of the softest photon while the softening begins at t c of the 
hardest photon; the ratio of the two times is the ratio of the photon 
energies reversed and squared, e.g., a factor of 1000 for the XRT 
band. 

The dependence of T a (E) on energy and grain size is 



T a (E) = r (E = lkcV,a = 0.1/im)(-^-)- s (- " 



-,4-9. 



;(5) 
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in the Rayleigh - Gans approximation, s — 2 (van de Hulst 
1957; Overbeck 1965; Hayakawa 1970; Alcock & Hatchett 1978; 
Mauche & Gorenstein 1986). 

The echo emission spectrum at an observer time t is obtained 

by 
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oc K ' I a 4 - q rt[x(E,a,t)]da. (6) 

The LC can be obtained by integrating Fe (t) over a desired energy 
bandpass. 

The softening can be seen from Eq. (6) as follows. Since 
x oc t}l 2 Ea, at some given time t, x might be > 1.5 for the 
hard photons and < 1.5 for the soft photons, while the interme- 
diate photon energy that defines and separates the "soft" and the 
"hard" corresponds to x ~ 1.5 and it decreases with time. For hard 
photons, j\ (x) oc x~ 2 . Taking E out of the integral in Eq. (6) gives 
Fs(t) oc S(E)E~ B ~ 2 , so the spectral index is increased by 4 for 
s = 2. For soft photons, j 2 (x) oc x 2 , so F E (t) oc S(E)E~ a+2 
and the spectral index is unchanged for s = 2. Therefore, the soft- 
ening happens first in the high energy part of the spectrum, and then 
propagates toward the lower energies with time, until the spectrum 
in the whole bandpass is softened - this is also when the plateau of 
the overall LC approaches to its end - with a change of the over- 
all spectral index A/3 = 4 with respect to the source spectrum (cf. 
Fig. 2 below). This change in /3 should be easy to detect if the echo 
emission dominates the plateau. 

We calculate the LCs of the dust echo in the XRT band for a 
variety of dust parameter values and the echo spectrum at different 
times. They are exactly the same as those obtained by Shao & Dai 
(2007) in their Fig. 3 and Fig. 4 and thus confirm their results and 
the analytical scalings derived above. 




t (s) 

Figure 1. The flux light curve in the XRT band of the echo emission from 
an extended dust zone with a distance range [R— , -R+]. The thick solid line 
(red in electronic version) is for the fiducial values used for the parameters 
of the dust spatial distribution model: R- = 50 pc, R+= 500 pc, 6= 2. The 
parameter values for the dust grain properties (same for all the curves) are 
a_ = 0.025 fim, a+ = 0.25 fim, q= 3.5, and s= 2. The parameter values listed 
in the legends are the only ones that are changed each time. The assumed 
source spectrum is a flat power law with a high energy cut-off: S(E) oc 
E° exp(— .E/100keV). For comparison, the LC from a single dust screen 
located at R = 100 pc with the same total optical depth is also shown here 
as the (black) thin solid line. 



2.1 Extended dust zone 

In the vicinity of GRBs, the dust zone may extend over a large 
distance. To study the difference in the echo emission properties 
of an extended dust zone and of a thin dust layer, we consider in 
this subsection a power-law dust distribution over a distance range 
[R-,R+] with the dust number density profile n(R) oc R~ s , 
where R is the distance to the source. The grain properties, e.g., 
size distribution, are assumed to be independent of R. For ease of 
calculation, we divide the extended dust zone into a series of N 
discrete thin dust layers (iV 30; a change of N does not affect 
the results), located progressively further from the GRB with equal 
separation in the log(R) scale. The scattered flux and its spectrum 
at any given time is the sum of the contributions from all dust layers 
at that time. 

Fig. 1 shows the LC of the scattered emission from an ex- 
tended dust zone for varied sets of parameters. It can be seen that 
the ending time of the plateau is mainly determined by the location 
of the inner boundary of the dust zone. This is not surprising be- 
cause the density of the dust is decreasing with radius thus the scat- 
tering LC arises mainly from the inner rim of the dust zone. Fig. 2 
shows the spectra of the scattered emission at different times, from 
which the softening is evident. The LCs in Fig. 1 and the spectra in 
Fig. 2 are almost same as the ones for a single dust layer model (cf. 
the thin solid line in Fig. 1 and Shao & Dai (2007)'s results in their 
Fig. 3 and 4), which shows that the generalisation of the model to 
an extended dust zone does not change much the temporal or spec- 
tral behavior of the scattered emission. 

We also calculate the overall spectral index in the XRT band 
at each observer time and the instantaneous decay index of the 
LC, which are shown in Fig. 3. Note that, due to the softening, 
the echo spectrum during the plateau is no longer a single power 
law function (see Fig. 2). Thus we calculate a "pseudo" spectral 
index /3n.3-io using the flux densities at the two ends of the XRT 
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Figure 2. The spectrum of the echo emission from an extended dust zone 
at different observer times. The propagation of the softening toward the low 
energies is evident. The model parameter values are the same as the fiducial 
ones in Fig. 1. The assumed source spectrum is a flat power law with a 
high energy cut-off: S(E) oc E° exp(-£/100keV). The change of the 
spectral index due to the softening is found to be insensitive to the source 
spectral index. 



band, 0.3 keV and 10 keV, respectively, to illustrate the extent of 
softening with respect to the source spectrum. 

The results show that the dust echo emission must experience 
significant spectral softening; the spectral slope increased by A/3 « 
3 from the early phase of the plateau to the end of the plateau. A 
more realistic extended dust zone model brings no notable change 
to this property. 
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Figure 3. The instantaneous flux decay index a (defined as 
i<0.3-iokevM ok t~ a ) and the spectral softening A/3 = 
A). 3- lOkcV ~~ Po f° r the echo emission from an extended dust zone, 
where /3 .3-iOkcV = log[F B (0.3keV)/F B (10keV)]/log(10/0.3) is 
the two-point spectral index and /3q is the source emission spectral index. 
In the cases plotted here we use /3rj = 0. By changing the values for /3g, 
e.g., to /3o= 1, -0.5 or -1, we find that the calculated A/3 is insensitive to 
/3rj. Each pair of lines for a and A/3 of the same line style correspond to a 
same set of model parameter values. The thick solid lines (red in electronic 
version) are for the same fiducial model parameter values used in Fig. 1. 
For comparison, the a and A/3 for the single dust screen model with the 
same parameter values as that in Fig. 1 are also shown as the thin solid 
lines (black in electronic version). 

3 SEARCH FOR SPECTRAL EVIDENCES IN THE DATA 

In this section we describe our search for the statistical evidence in 
the X-ray data during the plateau that can support the dust scatter- 
ing model. There are two pieces of evidence that we are looking for. 
First, if the plateau is due to the prompt X-rays scattering off the 
dust, in the early phase of the plateau when the spectral softening 
has not yet begun, the spectral index of the scattered emission must 
be the same as that of the prompt X-rays. Thus we expect to see 
in the data a correlation between the spectral index of the plateau, 
which we denote as p a here, and that of the prompt X-rays. There 
are two complications to note. (1) For the prompt emission, usually 
the X-ray spectral index is unavailable so we have to use the one 
for the prompt 7-rays, /3 7 , to represent it; in some cases the X-ray 
slope might be shallower than the 7-ray slope by 1 /2 due to a cool- 
ing break. (2) The published spectral index for the plateau is usually 
measured from the photon counts integrated over the whole plateau 
duration, therefore this spectrum might be softer than at the begin- 
ning of the plateau. Nevertheless, a mild trend of the correlation in 
the data should still be expected. 

The second evidence is based on the strong softening pre- 
dicted by the model as was demonstrated in Fig. 2-3, which show 
a strong evolution of the spectral index during the plateau and until 
its end. Thus, if the model is correct, the distributions of the spectral 
index during the plateau, p a , for a sample should be significantly 
smaller than that measured in the post-plateau phase, denoted as 

Pad- 

3.1 Sample 

A sample of GRBs showing X-ray plateaus with sufficient spec- 
tral and temporal information is needed to check for these two ev- 



idences. Willingale et al. (2007) analysed 107 Swift XRT detected 
GRB afterglows and found 80% of the bursts show a plateau in the 
X-ray LC. Out of the 80% of total bursts sample, 54 have both spec- 
tral indices before and after the end of plateau available. We further 
reduced the sample down to 26 bursts; we rejected those bursts that 
had one of the following properties: (1) the temporal decay slope 
a >0.8, too steep to be defined as a "plateau"; (2) XRT coverage is 
very sparse or long gaps exist during the plateau; (3) the "plateau" 
is actually due to one or more flares. We have also included 24 
bursts from the sample of Liang et al. (2007) that satisfy the above 
criteria. 

3.2 Results 

We plot /3 7 vs. p a and the difference between Pad and p a vs. p a 
in Fig. 4. No clear correlation between /3 7 and p a is seen, which 
disfavours the dust model. Generally p a is softer than /3 7 . This is 
consistent with the expected softening during the plateau. However, 
if it was the softening that could have weakened or broken the ex- 
pected correlation, there must be a bigger scatter in p a than in P-, . 
The data shows the contrary: p a is in the range of 0.5 - 1.5 where 
/3 7 is in -0.5 - 1.5. Thus the comparison between p a and /3 7 is in- 
consistent with the model expectation. 

Moreover, no dominant softening trend in the spectral index is 
seen from the plateau to the post-plateau phase; bursts with smaller 
Pa show slight softening and those with larger p a show slight hard- 
ening. Most of the bursts show zero spectral change across the 
end of the plateau within 1-<t measurement error. Only three bursts 
show evidence of spectral softening - two bursts at 2-a level (GRB 
050315: A/3= 0.4; GRB 060607A: Ap= 0.2) and one burst at 3-cr 
level (GRB 061202: A/?= 1.4). There are also two bursts showing 
spectral hardening - one at 2-cr level (GRB 060428A: A/3= -0.26) 
and one at 3-cr level (GRB 060413: Aj3= -0.93). Those individual 
cases are marked in the right panel of Fig. 4. 

The two results - no correlation between /3 7 and p a and no 
clear difference between p a and p a d, also reported in Willingale 
et al. (2007) and Liang et al. (2007) - are inconsistent with the 
expectations of the dust scattering model. 



4 TIME HISTORY OF THE SPECTRAL DURING THE 
PLATEAU 

The model predicts a significant spectral evolution from the begin- 
ning (~ 200 s after the burst) to the end of the dust echo plateau. 
The spectral index shows a monotonic increase by A/3 ~ 3 — 4. If 
the X-ray plateau is indeed due to or dominated by the dust echo 
emission, the strong spectral softening can be very easily detected 
in the XRT data. To compare the data with this expectation, we look 
closely at the time resolved spectral information during the plateau 
for the best observed GRBs. 

To determine whether there is any dominant trend of spec- 
tral evolution during the plateau phase, we compile a sample of 2 1 
GRBs with well defined plateau phases, excellent time coverage 
and good signal-to-noise ratios. For details of how the X-ray light 
curves used in this work were produced, see Evans et al. (2007). 
This sample is listed in Table 1. 

For each GRB, the overall XRT LC is considered to be com- 
posed of two components. The first one is the very rapid decay 
just following the 7-rays. The second component is the plateau and 
the subsequent normal power-law decay. Both components are well 
fitted by the same functional form as introduced by Willingale et 
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Shallow Phase Spectral Index ft Shallow Phase Spectral Index /S 

Figure 4. Left: The prompt phase BAT spectral index fi 1 vs. the shallow phase XRT spectral index a , Right: The change in the spectral index (3 aa — Pa 
vs. f3a, where /3 a( j is the spectral index of the decay after the end of the shallow phase. Several individual cases of GRBs which show evidence of spectral 
softening or hardening are labeled. The spectral index /3 is defined as /„ ex. v~ 13 . The error bars are 1-tr errors. The sample is selected from Willingale et al. 
(2007) (filled circles, and in black color in electronic version) and Liang et al. (2007) (filled triangles, and in red color in electronic version). 



al. (2007). We define Ti - the time of transition from the first to 
the second component - as the time when the two components are 
equal; this is a good measure of the start of the plateau. Tb is the 
end of the plateau and the start of the final power-law decay. 

We plot the hardness ratio, as defined by the ratio of the pho- 
ton counts in 1.5 - 10 keV and 0.3 - 1.5 keV bands, for each time 
interval of coverage during the plateau phase. We find that all the 
hardness ratio changes through the plateau phase are quite small 
and, for the bursts with the largest change the hardness ratios are 
getting harder (near the start of plateau), not softer. The ones that 
get softer do so only slightly. This confirms the findings by Butler 
& Kocevski (2007). 

We measure the spectral index taking into account photo- 
absorption at lower X-ray energies. We use the absorbed-power- 
law fit to the spectrum of the early XRT data - mainly the steep 
decline phase in the LC, which contains most of the photon counts 
- to determine the neutral H column density given as a combina- 
tion of two components - the Galactic column and a host intrinsic 
one. Then we use this neutral H absorption model to convert the 
measured hardness ratio in the plateau into the spectral index with 
an appropriate error. 

For many GRBs the time coverage during the plateau is rather 
patchy. Thus we select a time window which includes the plateau 
and takes into account the coverage. Sometimes it has to include 
data before Ti and after T2 so that the behaviour across the plateau 
is well constrained. The evolution of the spectral index over the 
selected time window is fitted as a linear function of the logarithmic 
time. Extrapolating the best fit function to both sides of the window 
gives the spectral indices /3i and 02 at Ti and T2, respectively. 
These are the best estimates of the spectral index at the start and 
the end of the plateau. 

Fig. 5 shows a few examples of the /9-evolution during the 
plateau. The observed /9-evolutions for our sample are tabulated 
in Tab. 1. None of the afterglows show a notable softening over 
the plateau. Most show zero evolution of f3, with very small uncer- 
tainty. A few show a small, marginally significant, hardening. For 
the examples shown in Fig. 5, we also add in the lower panels of 
Fig. 5 the expected /3-evolutions from the dust scattering model. 
A few model parameters (R, q and s) were set free to change and 



then were optimised in each example in order to best reproduce the 
observed plateau and post-plateau LC. The expected is system- 
atically larger than the observed one even at the beginning of the 
plateau phase, because the softening has already begun there. The 
expected strong evolution in j3 distinctly differs from the stableness 
of the observed /3 during the plateau. 

To summarise this section, though the dust scattering model 
can nicely fit the LCs of plateau and post-plateau decay (see also 
Shao et al. 2008), the expected large value and strong evolution of 
the spectral index sharply contradict the data. 



5 X-RAY LIGHT CURVES IN SOFT AND HARD 
ENERGY CHANNELS 

Another prediction from the dust scattering model is the different 
temporal behaviours of the LCs in low and high energy channels. 
This can be seen from equation (4): observations carried out at pho- 
ton energy E should see the end of the plateau at t c oc E~ 2 . This 
feature is demonstrated in the top panel of Fig. 6 where we calcu- 
lated the LCs of scattered emission from an extended dust zone in 
0.3 - 1.5 keV and 1.5 - 10 keV, respectively. It shows that the soft 
X-ray LC has a more extended plateau than the hard X-ray one. 

We find that for all bursts in our sample of plateau GRBs, the 
temporal behaviours in soft vs. hard X-ray LCs look identical. As 
an example, the bottom panel of Fig. 6 shows the soft vs. hard X- 
ray LCs of GRB 061 121 which has a long, dense time coverage and 
the best photon statistics among all bursts in the sample (Page et al. 
2007). The soft and hard X-ray LCs for this burst are identical. This 
feature rules out the dust scattering model for the plateau. 



6 OPTICAL EXTINCTION DUE TO THE DUST 

The dust grains which scatter the X-ray photons will also cause 
extinction in the optical band. This can provide an additional con- 
straint for the dust scattering model. Thus we estimate the extinc- 
tion in V Band, Ay, caused by the dust required for the model and 
compare it with Ay derived directly from optical observations. 
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Figure 5. The X-ray LC and the time history of the spectral index j3 for four examples of GRBs with plateaus and the corresponding expectations from the 
dust scattering model. In the upper panel for each example, Ti and T2 mark the beginning and the end of the plateau phase, respectively. The plateau phase 
and the post-plateau decay LC is mimicked by the dust scattered emission through adjusting the parameters (R, q and s) of a single dust screen model (using 
an extended dust zone model does not change the result), shown as the solid line. In the lower panel for each example, the time-resolved /3 are plotted as filled 
circles inside the [Ti, T2] window and as crosses outside the window. A linear function of log(t), as shown by the solid straight line, is fitted to the time 
history of /3 within the [Ti, T2] window. The fit extrapolation at Ti and T2 gives f3\ and /?2, respectively, as marked by the filled bars. The dashed lines are 
the expected /3-evolutions for the dust models that were optimised in the upper panels. 



Predehl & Schmitt (1995) found an empirical relation between 
the X-ray dust scattering optical depth t(E) and Av for the X-ray 
halos of 24 galactic X-ray point sources: 

t{E) = 0.06A v (£/lkeV)" 2 , (7) 

where t(E) is obtained from modelling the X-ray halo surface 
brightness distributions with dust grain properties similar to the 
fiducial ones we used in our calculations. Draine & Bond (2004) 
derived a similar relation based on a dust model developed by 
Draine (2003): 

t(E) = 0.15Av(£/lkeV) -1 - 8 . (8) 

Note that Eq. (7) is for the dust in the Milky Way (MW), while for 
GRB hosts most absorption fits tend to favour the Small Magellanic 
Cloud (SMC) extinction law (e.g., Schady et al. 2007; but see the 
discussion below toward the end of this section). It was shown that 
the difference between the the MW and SMC extinction laws can be 
well reproduced in a model by adjusting the relative abundances of 
graphite and silicate grains, while leaving all other dust properties 
fixed; in this case Av and t(E) at E ^ 7 eV are both the same for 
these two environments (Pei 1992). Thus it is viable to apply Eq. 
(7) to GRB hosts. 

The above relations are for quantities in the rest frame of the 
source. If the source is at cosmological distances, like GRB hosts, 
we have to take into account the cosmological redshift of the photon 
energy when calculating t(E) and Av from the observed quanti- 



ties. The dust scattering optical depth, to = r(E — IkeV), in the 
rest frame of the GRB host at a redshift 2 can be estimated by 

Sx.i ~ (l + 2)" 2 T S 7 ,i, (9) 

where S 7 ,i = S 1 (E = IkeV) is the specific fluence extrapolated 
from the Burst Alert Telescope (BAT) total fluence during the burst 
to 1 keV, Sx,i = Sx(E — IkeV) is the specific fluence at 1 keV 
during the plateau phase, both measured in the observer's frame. 
There is a factor of (1 + z)~ 2 because the ratio (Sx,i/Sj,i) is 
actually equal to the host rest-frame t(E) at E = (1 + z) keV and 
t(E) <x E~ 2 (cf. Eq. 5). 

We select a sub-sample of GRBs which have good XRT tem- 
poral coverage during the afterglow phase from the sample of Liang 
et al. (2007) that provided S-, and Sx for those bursts. Then to is 
calculated and the associated Av is inferred from to via Eq. (7) 
and (8). The sample and the results are listed in Table 2. Almost in 
all cases (except for two) to is > 1 and some even have to > 10, 
which means this model requires that only a very tiny fraction of 
photons with energy of 1 keV in the host rest frame can escape the 
dust without scattering. None of the sub-sample have Av < 1 and 
85% of them have A v > 10. 

In comparison, Schady et al. (2007) determined Av for 6 GRB 
afterglows from the Ultraviolet and Optical Telescope (UVOT) to 
XRT Spectral Energy Distributions (SEDs), with Av ranging from 
0.1 to 0.7. These extinctions are significantly smaller than that ex- 
pected from the dust scattering model. Note that a considerable 
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Table 1. The spectral indices at the beginning and at the end of the X-ray plateaus for a sample of GRBs. From a larger sample of GRBs showing well defined 
plateau phases, only those with high signal-to-noise ratio and long time coverage in the plateau are selected. T\ marks the transition of the LC from the prompt 
component to the plateau component. T2 marks the end of the plateau and the transition to the final power-law decay. @\ and ft are given by the extrapolation 
of a function fit to the evolution of the spectral index of available data within a window defined by T\ and T2 . The error in [3 is at 90% confidence level. 
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Table 2. The fiducial host rest frame optical depth at E = 1 keV required by the dust scattering model and the associated rest frame visual extinction for a 
sub-sample of GRBs with well observed shallow X-ray decays. The fluence and spectral data are from Liang et al. (2007). Ay is given via equation (7) and 

via equation (8). References to the GRB redshifts: 050315 - Kelson & Berger (2005); 050319 - lakobsson et al. (2006); 050401 - Fynbo et al. (2005); 
050803 - Bloom et al. (2005); 060210 - Cucchiara et al. (2006); 060714 - lakobsson et al. (2006); 060729 - Thoene et al. (2006); 060814 - Thoene et al. 
(2007). z = 2 is assumed for GRBs without known z. 
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fraction (1/4 - 1/3) of the plateau X-ray afterglows have bright op- 
tical counterparts (Figure 2 of Liang et al. 2007). 

But a cautiousness has to be taken regarding the Av deter- 
minations above. There are some recent studies that show the dust 
properties of GRB hosts does not resemble those of any galaxy in 
our neighborhood. In particular, Chen et al. (2006), Perley et al. 
(2008) and Li et al. (2008) have found for some GRBs the mod- 
eled extinction curve is "gray", i. e., much flatter than any of the 
templates (MW, SMC, the Large Magellanic Cloud). Since our de- 
termination of Av has used the Av - r relations appropriate for 



these template-type dust, these findings are likely to raise uncer- 
tainties in the determined Av- 

The dust grain size distribution is usually described by a power 
law with an index q and within a range of grain size (a_, a + ). A 
"gray" extinction curve could be due to a flatter grain size distribu- 
tion (smaller q) or a larger a+, as suggested by Li et al. (2008). We 
calculated Av and to independently for a dust model with the com- 
position resembling the SMC and the "gray" type that these authors 
have found for some GRBs, to see the dependence of their ratio on 
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Figure 6. The LCs in the soft and the hard X-ray channels, respectively, 
as predicted by the dust echo model and observed in GRBs. Top: The LCs 
of the echo emission from an extended dust zone model. The values of the 
model parameters R— and are chosen such that the echo emission 0.3 
- 1 keV LC best mimics the plateau of GRB06 1121. Other model parame- 
ters are the same as the fiducial ones in Fig. 1 . Bottom: The LCs of a typical 
GRB 'plateau' in 0.3 - 1.5 keV (squares) and 1.5 - 10 keV (triangles), re- 
spectively. The early rapid decay at t < 200 s of the observed LC does not 
belong to the plateau phase and it is thought to be of a different origin. 



the grain size distribution parameters. The Av is calculated by the 
following equation (Weingartner & Draine 2001) 



Av = (2.57rloge) 



ex t{a,\v)^r^-da, (10) 
da 



where Q ex t is the extinction efficiency factor, usually a function 
of grain size and photon energy. Draine & Lee (1984) calculated 
Q ex t(a,\) for graphite and silicate grains. Since the GRB dust 
environments are described either by SMC or "gray" type extinc- 
tion curve (Shady et al. 2007; Li et al. 2008), for both of which a 
good dust composition model needs silicate only (Pei 1992; Li et 
al. 2008), we use Q e xt(a, A) for silicate only. The dust scattering 
optical depth can be calculated by to = r a (E = lkeV)da, 
where T a (E) is given by Eq. (5). Note that the normalizations of 
Ay and to are unimportant here because we are looking at only 
their dependences on a+ and q. 

The ratio of Av/tq for varying a+ and q is shown in Fig. 
7. We find that Av/tq is only slightly dependent on q - Av/to 
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Figure 7. The dependence of the ratio Av/tq of the dust on the grain size 
distribution parameters a+ and q. The dust is assumed to have a grain size 
distribution dN(a)/da on a~ q within the grain size range [a_ , o+], where 
AT (a) is the column density of all grains with size a. 



decreases by factor of ~ 2 for q changing from 3.5 to 2.6; but 
Av/to is more sensitive to a+ - it decreases by a factor of ~ 10 
when a+ changes from 0.25 fim to 2.5 /im. Even after taking these 
effects into account, the Av we obtained for our GRB sample are 
still very large (for one half of the sample, Av > 10). One of the 
real problems with the dust scattering model for the X-ray plateau 
phase is that it requires to ^ 10 for one half of our sample (see 
Table 2), which seems physically unreasonable. 

Thus, the dust scattering model is not a viable explanation for 
the X-ray plateaus because of the large extinction in the optical 
band it predicts but not observed. 



7 CONCLUSION AND DISCUSSION 

We have shown that in the dust scattering model the scattered X- 
ray emission must experience strong softening spectral evolution, 
with a significant change of the spectral index in 0.3 - 10 keV of 
A/3 ~ 2 — 3 from the emerging of the plateau to its end. However, 
for a sample of GRBs with X-ray plateaus and with good quality 
data, no softening spectral evolution during the plateau phase is 
found, and in a few cases even traces of slight hardening are seen. 

The change of (5 according to the model does not depend on 
the spectral index of the source emission. The Rayleigh - Gans ap- 
proximation is used in this paper to calculate the scattering cross 
section of the dust grain. It was claimed that this approximation 
tends to overestimate the scattering efficiency below 1 keV, typi- 
cally by a factor of 4 at 0.5 keV and a factor of 2 at 1 keV, mainly 
due to the absorption of the soft X-ray photons by the K and L shell 
electrons in the dust grain (Smith & Dwek 1998), and that could 
change the spectral slope at the soft end (< 1 keV) and counteract 
against the softening (Shao et al. 2008). But we argue that this ef- 
fect dose not alleviate the expected softening, because the discrep- 
ancy between the real scattering cross section and the Rayleigh - 
Gans approximation caused by this effect, which is mainly below 1 
keV, must have been largely accounted for by the required neutral 
H absorption in the routine power-law fit to the plateau spectra. 
The XRT spectral index is mainly determined by the photons with 
energy above 1 keV which is not affected by this effect. Moreover, 
this effect is time independent while the softening we consider is 
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a strongly time dependent behaviour. Dust destruction by the GRB 
prompt emission is of very little relevance here because it happens 
within a distance smaller than the location of the dust considered 
in this work (e.g., Waxman & Draine 2000). Thus the Rayleigh - 
Gans approximation is sufficiently accurate for the effect consid- 
ered in this work. 

The dust scattering model also predicts very different tempo- 
ral behaviours in the soft X-ray vs. hard X-ray LCs; the plateau 
lasts longer in soft X-rays. But this feature is not found in the data. 
Furthermore, the large scattering optical depth of the dust required 
by this model in order to explain the X-ray plateaus leads to ex- 
tremely large extinction in optical - Av > 10. This is inconsistent 
with the observed extinctions for GRBs. 

We conclude that the dust scattering model, though very at- 
tractive, can not explain the X-ray plateaus seen in most GRB af- 
terglows. Although it is very likely that dust exists near the site of 
GRBs, and will scatter some fraction of the prompt and afterglow 
X-rays, this scattered emission is not a dominant contributor to the 
observed X-ray plateau. For those cases where an achromatic break 
at the end of the plateau is seen, a late, steady energy injection to 
the external shock is a more likely mechanism for producing the 
observed X-ray plateau, though it may not be able to work well for 
the cases with chromatic breaks. 
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